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1. Introduction 


High energy experiments like Tevatron or LHC study the Standard Model of particle physics 
and its possible extensions. At the core of these experimental studies are Monte Carlo (MC) gen¬ 
erators that are based on factorization theorems. Their construction usually involves additional 
approximations and phenomenological models. However, this approach allows the simulation of 
particle scatterings starting from the har'd interaction and then dressing the external legs with further 
radiation generated by the parton shower. This will evolve the hard state down to a low energy scale 
where non-perturbative effects are important. At that point hadronization models are employed. 

In order to keep up with the increasing precision of experimental data, MC generators have to 
be improved. There are several options to achieve this. For fixed-order calculations the inclusion 
of next-to-leading order (NLO) corrections in quantum chromodynamics (QCD) are widely auto¬ 
mated ||T} Currently, efforts are ongoing in the automatization of electroweak corrections [|1] 
while NNLO QCD corrections are only available for 2 —)• 2 processes. 

Fixed-order calculations usually suffer from large logarithms, which can be resummed to all 
orders in perturbation theory using analytical resummation 0 or parton shower methods liiSi- 
Matching of NLO fixed-order calculations fo parton showers has the benefit of combining several 
features. Well separated partons are then correctly described by matrix elements in perturbative 
quantum field theory, whereas the soft and collinear parton splittings are generated by parton show¬ 
ers. There are several methods in the literature to do the matching, the most popular ones being 
Mc@Nlo and Powheg [^. The current matching and merging schemes are all limited by the 
accuracy of the shower algorithms, that include only leading colour (LC) and leading logarithmic 
(LL) accuracies and no spin correlations. In order to go beyond these approximations one has to 


include soft-gluon interferences and other subleading effects. In Ref. |10] it has been shown that 
these subleading effects can be sizable for specific observables. 

In this proceeding, we first briefly review fhe basics of fhe Nagy-Soper shower. In secfion 3 
we explain fhe mafching scheme for fhis shower and in secfion 4 we show first results for the pp ^ 
ttj + X production at the LHC. Finally, we conclude and give an outlook for future improvements 
of this work in section 5. 


2. Nagy-Soper parton shower 


Here, we give a very brief summary of the Nagy-Soper parton shower introduced by Zoltan 
Nagy and Davison Soper 00- We only highlight the necessary concepts to understand the 
parton shower matching, while a more thorough discussion can be found in Ref. [[1^]. We start 
from the all-order expression of the expectation value of an observable F, for a 2 —m process 


111 J 


where fi^Sf^Ci and pi represent the flavour, spin, colour and momentum of a particle. The general- 
ized phase space integration measure, [d{p,f,s',c',s,c},n], includes the integration over the initial 
state momentum fractions rja and rjh as well as the summation over spin and colour indices. The 
matrix element, |-#({p,/}m))^ is a vector in colour and spin space and its square can be rewritten 
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in the form of a quantum density matrix defined as 




fa{Pa,pj)fb{Pb,Pj) 


Anc{a)nc{b) x flux 

( 2 . 2 ) 

The parton shower evolves this quantum density from the hard scale, to, of the scattering process 
down to the low scale, tf, where non-perturbative physics models can be applied. The evolution of 
the system is described by a unitary operator U {tf,to)-> which satisfies 


dU [tfo) 
dt 


= [jif,{t)-'r{t)]u{t,to) 


(2.3) 


Here fhe evolution operator U(tffo) is given by the real splitting operator, Jifjf), which describes 
the emission of resolved particles and the virtual operator, Y(t), which describes unresolved emis¬ 
sions. The above differential equation is solved by the following operator 


ld{tF,to) =N{tf,to)+ [ dT:[Jfi{x)— 'fs{'^)]U{T:,to) 

Jtn 


with the normal Sudakov form factor 


= exp r/Ti^(T)^ . 


(2.4) 


(2.5) 


The virtual splitting operator is decomposed into "^(t) = + 'i^(t), where "^(f) is the colour 

diagonal part and the colour off-diagonal part with subleading contributions. The colour 
diagonal part will be exponentiated while the off-diagonal contribution is treated as perturbation. 

This shower concept differs in many ways from conventional showers and has been partly 
implemented in a MC program Deductor [[I^. A main conceptual difference is that the splitting 
functions are derived using factorization on the amplitude level [|ri|]. In addition, initial state charm 
and bottom quarks are treated as massive and PDFs used in the shower are evolved according to 
the shower splitting kernels [^, |T^]. The shower is consistently able to include spin and colour 
correlations throughout the whole parton evolution [12, |I^, [l^. The ordering parameter in the 
evolution is a virtuality based one and implements the validity of the on-shell approximation in 


each step of the evolution [18] 


e ’ = — 


^2 ^ |(p/±Pm+l) -»r(0| g2 


( 2 . 6 ) 


2prQ 

Finally, the Nagy-Soper shower uses a global momentum mapping that improves resummation 


effects in Drell-Yan Z-production [19]. 


3. Parton shower matching 


At next-to-leading order the quantum density matrix has to be extended by including the virtual 
and real matrix elements 


Ip)— |pi^ + + Ip^+i) 

Born, ^(1) Virtual, Real, 


,(l)i 


,(0) 


(3.1) 
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Applying the shower evolution operator to this density generates a spurious non-zero contribution 
at which can be easily seen by expanding the evolution operator 

iF\UitF,to)\p)^iF\p)+ r dr {F\[J^,{r)-yir)]\pj^^) +. (3.2) 

Jta 

As a consequence, the total cross section is changed and the first emission is double counted. In 
order to remove this additional contribution the Mc@Nlo approach has been applied. We will first 
focus on fully inclusive processes like pp ^tt and then we will explain the matching of exclusive 
processes which already suffer from divergences at the LO. 


3.1 Fully inclusive processes 

The idea is to redefine fhe quanfum densify mafrix by providing suifable counferferms for fhe 
shower confribufion 

Ip) = 1p)- f " dT [ j ^ f ’ j { T )-' r { T )]\ p ! n '^) + ff { aj ) . (3.3) 

Jto 


By dropping fhe infrared cufoff {tf —)• oo) we observe fhaf fhe shower nafurally incorporafes fhe 
NLO subfracfion scheme 


noo poo 

/ r/TJ^(T) = £S// r/T5(T-h)0(T-to)=£S/0(h-to), 

Jt(\ 1 Jtp 1 



-Viz) = £ / dr;S/0(h - to) = I(to) + K(to) 


(3.4) 


Thus, we see fhaf mafching is a fwo-sfep procedure 


a[F] = -^1 [d<I>,„](f |[/(tf,to)|<I>,„)(<I>m|5) 

+ . ] [[d^m+l]{F\UitF,to)\^m+i)i^m+l\H) , 

(ni+ 1)! J 

where d>,„ = {p,f,c,s,c,s}m- Firsf we have fo generafe fhe samples 


(3.5) 


(d>„,|S) = (<I>,„|pi°)) + i%n\pln^) + (<!>„, |I(to) +K(to) +P|p,lf^) , 
i^,n+i\H) = (<!>„,+i|p,S2i) -£(d>„,+i|S,|p,lf^)0(t/ -fo) , 

I 

and fhen apply fhe shower evolufion operator U[tF,to)- 


3.2 Exclusive processes 

The sifuafion is more complex for exclusive processes. In order fo avoid double counfing, we 
have fo add inclusive jef functions Fj and we have fo modify fhe subfracfion ferms 


a[F] = ^J [d<^nr]iF\UitF,to)\<^,n)i<^,n\S)Fii{pJU) 

+ [f/‘I>„,+i](F|f/(ff,fo)|<I>™+i)(4>,„+i|.y)F/({p,/Wi) , 


(3.7) 
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with 

(<I>„+i|^) = -£(4>,„+i|S/|pif))0(h -to)i7(<2/({p , (3.8) 

where the jet function F[{Qi{{p,f},„^i)) = Fi{{p,f},„) acts on the underlying Born kinematics. 
Expanding the evolution operator now gives 

.] m\ {m+l)\JtQ ( 39 ) 

X (d>„,|pi°^) [1 -F7({p,/Wi)]F;({p,/},„) + . 

Thus, double counting is removed if E/({p,/},„+i) = 1 for E({p,/}m+i) / 0, i.e. when generation 
cuts are more inclusive than cuts on the final observable. 


4. Application: pp ^ ttj +X 


The implementation details of the scheme presented in the previous section in the Helac- 
Nlo framework [ |20| ] can be found in Ref. [Oi- We now present results for pp —)• tlj + X at the 
LHC. The NLO QCD corrections have been already presented in Ref. [^]. Lirst parton shower 
matched calculations using the POWHEG method were presented in Ref. |^], while merging sev¬ 
eral matched calculations for different jet multiplicities is discussed in Ref. [23]. 

We consider the LHC at y/s = 8 TeV. The top quark mass is m, = 173.5 GeV. In the shower we 
set the parton masses of charm and bottom quarks to nic = 1.4 GeV and nth = 4.75 GeV. We use the 
MSTW2008NLO PDL [24] set in our calculation and provided it at pf = I GeV to Deductor. 
Renormalization and factorization scales are set to pn = pf = po = nif. We use the anti-kr jet 
algorithm [^] with AR = 1 and the analysis cuts are ptUi) > 50 GeV and |y(7i)| < 5, while the 
generation cut is pr( 7 i) > 30 GeV. The initial shower time is chosen to be 


e ‘° = min 

¥j 


'^PrPj\ 


(4.1) 


where the parameter pr allows us to address parton shower uncertainties. Lor our central prediction 
we choose pr = 1. The shower evolution is restricted to LC and spin-averaged contributions. In 
addition we do not include non-perturbative effects or top quark decays. The comparison with other 
Monte Carlo generators is, therefore, performed at the level of the perturbative evolution. Lor the 
comparison we use aMc@NLO [Q] with the Mc@Nlo matching in conjunction with Pythia8 
[§] (pr ordered parton shower) and Pythia6Q [Q] (virtuality ordered parton shower). Lor the 
PowHEG matching we use the PowHEG-Box [ ^^ results together with Pythia8. 

In Pig. H] we address scale and parton shower uncertainties: /2 < po < Inif and 1/2 < < 

2. The pt of the hardest jet shows a flat and reduced scale variation with respect to both parameters 
Po and Pt- This is expected since it is already an NLO-accurate observable. On the contrary, the 
Pt of the ttj\ system presents a stronger dependence on the shower parameter pr in the low pj 
regime, where the parton shower dominates. The scale dependence sets in once the real matrix 
element is present, i.e roughly at 30 GeV, and grows rapidly with transverse momentum. 

Pig. ^ illustrates the comparison with other MC generators. Lor inclusive observables, like 
the Pt of the top quark in the left panel, we do not observe substantial differences between the 
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MC generators. This is expected because this observable is already accurate at NLO and serves 
as a good cross check of our implementation. For exclusive observables like the pj of the ttj 
system, that is shown in the right panel of Fig. ^ we find a strong dependence on the initial shower 
conditions. We observe that aMc@NLO+PYTHlA and Powheg+Pythia overshoot the high px 
tail, where one would like to recover the real matrix element description. On the other hand, 
Helac-Nlo+Deductor and aMc@NLO+PYTHlA6Q preserve the prediction in the high energy 
tail. The comparison for other observables can be found in Ref. [13]. 



Figure 1: Differential cross section as a function of the transverse momentum of the hardest jet (left panel) 
and of the ttji system (right panel). 




Figure 2: Differential cross section as a function of the transverse momentum of the top quark (left panel) 
and of transverse momentum of the ttji system (right panel). 
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5. Summary and Outlook 

We presented the NLO matching scheme for the Nagy-Soper parton shower in the spirit of the 
Mc@Nlo method. We studied the production of tlj at the LHC using the Helac-Nlo-i-Deductor 
framework and compared it to the other MC generators. We want to stress that the current accuracy 
of Helac-Nlo-i-Deductor is only LC and spin-averaged, i.e the same as for other MC pro¬ 
grams. However, this comparison is an important validation of our implementation. Differences, 
present in exclusive observables, can be traced back, for example, to the choice of the initial shower 
conditions. 

In the future we want to extend the matching scheme implemented in Helac-Nlo to include 
full colour and spin correlations. Nevertheless, those effects have to be first added to Deductor. 
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